Distinct structural variants and repeat landscape shape the genomes of the ancient grapes Aglianico and Falanghina.

Mounting evidence recognizes structural variations (SVs) and repetitive DNA sequences as crucial players in shaping the existing grape phenotypic diversity at intra- and inter-species levels. To deepen our understanding on the abundance, diversity, and distribution of SVs and repetitive DNAs, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), we re-sequenced the genomes of the ancient grapes Aglianico and Falanghina. The analysis of large copy number variants (CNVs) detected candidate polymorphic genes that are involved in the enological features of these varieties. In a comparative analysis of Aglianico and Falanghina sequences with 21 publicly available genomes of cultivated grapes, we provided a genome-wide annotation of grape TEs at the lineage level. We disclosed that at least two main clusters of grape cultivars could be identified based on the TEs content. Multiple TEs families appeared either significantly enriched or depleted. In addition, in silico and cytological analyses provided evidence for a diverse chromosomal distribution of several satellite repeats between Aglianico, Falanghina, and other grapes. Overall, our data further improved our understanding of the intricate grape diversity held by two Italian traditional varieties, unveiling a pool of unique candidate genes never so far exploited in breeding for improved fruit quality.

Structural Variation Evolution at the 15q11-q13 Disease-Associated Locus.

The impact of segmental duplications on human evolution and disease is only just starting to unfold, thanks to advancements in sequencing technologies that allow for their discovery and precise genotyping. The 15q11-q13 locus is a hotspot of recurrent copy number variation associated with Prader-Willi/Angelman syndromes, developmental delay, autism, and epilepsy and is mediated by complex segmental duplications, many of which arose recently during evolution. To gain insight into the instability of this region, we characterized its architecture in human and nonhuman primates, reconstructing the evolutionary history of five different inversions that rearranged the region in different species primarily by accumulation of segmental duplications. Comparative analysis of human and nonhuman primate duplication structures suggests a human-specific gain of directly oriented duplications in the regions flanking the GOLGA cores and HERC segmental duplications, representing potential genomic drivers for the human-specific expansions. The increasing complexity of segmental duplication organization over the course of evolution underlies its association with human susceptibility to recurrent disease-associated rearrangements.

Occurrence of L1M Elements in Chromosomal Rearrangements Associated to Chronic Myeloid Leukemia (CML): Insights from Patient-Specific Breakpoints Characterization.

Chronic myeloid leukemia (CML) is a rare myeloproliferative disorder caused by the reciprocal translocation t(9;22)(q34;q11) in hematopoietic stem cells (HSCs). This chromosomal translocation results in the formation of an extra-short chromosome 22, called a Philadelphia chromosome (Ph), containing the BCR-ABL1 fusion gene responsible for the expression of a constitutively active tyrosine kinase that causes uncontrolled growth and replication of leukemic cells. Mechanisms behind the formation of this chromosomal rearrangement are not well known, even if, as observed in tumors, repetitive DNA may be involved as core elements in chromosomal rearrangements. We have participated in the explorative investigations of the PhilosoPhi34 study to evaluate residual Ph+ cells in patients with negative FISH analysis on CD34+/lin- cells with gDNA qPCR. Using targeted next-generation deep sequencing strategies, we analyzed the genomic region around the t(9;22) translocations of 82 CML patients and one CML cell line and assessed the relevance of interspersed repeat elements at breakpoints (BP). We found a statistically higher presence of LINE elements, in particular belonging to the subfamily L1M, in BP cluster regions of both chromosome 22 and 9 compared to the whole human genome. These data suggest that L1M elements could be potential drivers of t(9;22) translocation leading to the generation of the BCR-ABL1 chimeric gene and the expression of the active BCR-ABL1-controlled tyrosine kinase chimeric protein responsible for CML.

A high-resolution map of small-scale inversions in the gibbon genome.

Gibbons are the most speciose family of living apes, characterized by a diverse chromosome number and rapid rate of large-scale rearrangements. Here we performed single-cell template strand sequencing (Strand-seq), molecular cytogenetics, and deep in silico analysis of a southern white-cheeked gibbon genome, providing the first comprehensive map of 238 previously hidden small-scale inversions. We determined that more than half are gibbon specific, at least fivefold higher than shown for other primate lineage-specific inversions, with a significantly high number of small heterozygous inversions, suggesting that accelerated evolution of inversions may have played a role in the high sympatric diversity of gibbons. Although the precise mechanisms underlying these inversions are not yet understood, it is clear that segmental duplication-mediated NAHR only accounts for a small fraction of events. Several genomic features, including gene density and repeat (e.g., LINE-1) content, might render these regions more break-prone and susceptible to inversion formation. In the attempt to characterize interspecific variation between southern and northern white-cheeked gibbons, we identify several large assembly errors in the current GGSC Nleu3.0/nomLeu3 reference genome comprising more than 49 megabases of DNA. Finally, we provide a list of 182 candidate genes potentially involved in gibbon diversification and speciation.

Yeast as a Model to Unravel New BRCA2 Functions in Cell Metabolism.

Mutations in BRCA2 gene increase the risk for breast cancer and for other cancer types, including pancreatic and prostate cancer. Since its first identification as an oncosupressor in 1995, the best-characterized function of BRCA2 is in the repair of DNA double-strand breaks (DSBs) by homologous recombination. BRCA2 directly interacts with both RAD51 and single-stranded DNA, mediating loading of RAD51 recombinase to sites of single-stranded DNA. In the absence of an efficient homologous recombination pathway, DSBs accumulate resulting in genome instability, thus supporting tumorigenesis. Yet the precise mechanism by which BRCA2 exerts its tumor suppressor function remains unclear. BRCA2 has also been involved in other biological functions including protection of telomere integrity and stalled replication forks, cell cycle progression, transcriptional control and mitophagy. Recently, we and others have reported a role of BRCA2 in modulating cell death programs through a molecular mechanism conserved in yeast and mammals. Here we hypothesize that BRCA2 is a multifunctional protein which exerts specific functions depending on cell stress response pathway. Based on a differential RNA sequencing analysis carried out on yeast cells either growing or undergoing a regulated cell death process, either in the absence or in the presence of BRCA2, we suggest that BRCA2 causes central carbon metabolism reprogramming in response to death stimuli and encourage further investigation on the role of metabolic reprogramming in BRCA2 oncosuppressive function.

BL1391: an established cell line from a human malignant peripheral nerve sheath tumor with unique genomic features.

Malignant peripheral nerve sheath tumors (MPNSTs) are aggressive tumors, accounting for around 5% of all soft tissue sarcomas. A better understanding of the pathogenesis of these tumors and the development of effective treatments are needed. In this context, established tumor cell lines can be very informative, as they may be used for in-depth molecular analyses and improvement of treatment strategies. Here, we present the genomic and transcriptomic profiling analysis of a MPNST cell line (BL1391) that was spontaneously established in our laboratory from a primary tumor that had not been exposed to genotoxic treatment. This cell line shows peculiar genetic features, such as a large marker chromosome composed of high-copy number amplifications of regions from chromosomes 1 and 11 with an embedded neocentromere. Moreover, the transcriptome profiling revealed the presence of several fusion transcripts involving the CACHD1, TNMA4, MDM4, and YAP1 genes, all of which map to the amplified regions of the marker. BL1391 could be a useful tool to study genomic amplifications and neocentromere seeding in MPNSTs and to develop new therapeutic strategies.

Evolutionary Dynamics of the POTE Gene Family in Human and Nonhuman Primates.

POTE (prostate, ovary, testis, and placenta expressed) genes belong to a primate-specific gene family expressed in prostate, ovary, and testis as well as in several cancers including breast, prostate, and lung cancers. Due to their tumor-specific expression, POTEs are potential oncogenes, therapeutic targets, and biomarkers for these malignancies. This gene family maps within human and primate segmental duplications with a copy number ranging from two to 14 in different species. Due to the high sequence identity among the gene copies, specific efforts are needed to assemble these loci in order to correctly define the organization and evolution of the gene family. Using single-molecule, real-time (SMRT) sequencing, in silico analyses, and molecular cytogenetics, we characterized the structure, copy number, and chromosomal distribution of the POTE genes, as well as their expression in normal and disease tissues, and provided a comparative analysis of the POTE organization and gene structure in primate genomes. We were able, for the first time, to de novo sequence and assemble a POTE tandem duplication in marmoset that is misassembled and collapsed in the reference genome, thus revealing the presence of a second POTE copy. Taken together, our findings provide comprehensive insights into the evolutionary dynamics of the primate-specific POTE gene family, involving gene duplications, deletions, and long interspersed nuclear element (LINE) transpositions to explain the actual repertoire of these genes in human and primate genomes.

1q23.1 homozygous deletion and downregulation of Fc receptor-like family genes confer poor prognosis in chronic lymphocytic leukemia.

The identification of chromosome 1 translocations and deletions is a rare and poorly investigated event in chronic lymphocytic leukemia (CLL). Nevertheless, the identification of novel additional molecular alterations is of great interest, opening to new prognostic and therapeutic strategies for such heterogeneous hematological disease. We here describe a patient affected by CLL with a mutated IGHV status, showing a balanced t(1;3)(q23.1;q21.3) translocation and a der(18)t(1;18)(q24.2;p11.32), accompanying the recurrent 13q14 heterozygous deletion in all analyzed cells at onset. By combining whole-genome sequencing, SNP array, RNA sequencing, and FISH analyses, we defined a 1q23.1 biallelic minimally deleted region flanking translocations breakpoints at both derivative chromosome 1 homologues. The deletion resulted in the downregulation of the Fc receptor-like family genes FCRL1, FCRL2, and FCRL3 and in the lack of expression of FCRL5, observed by RT-qPCR. The mutational status of TP53, NOTCH1, SF3B1, MYD88, FBXW7, and XPO1 was investigated by targeted next-generation sequencing, detecting a frameshift deletion within NOTCH1 (c.7544_7545delCT). We hypothesize a loss of tumor suppressor function for FCRL genes, cooperating with NOTCH1 mutation and 13q14 genomic loss in our patient, both conferring a negative prognosis, independently from the known biological prognostic factors of CLL.

Correction: MYC-containing amplicons in acute myeloid leukemia: genomic structures, evolution, and transcriptional consequences.

In the original version of this Article, the affiliation details for Giovanni Martinelli were incorrectly given as 'Department of Experimental, Diagnostic and Specialty Medicine, University of Bologna, 40138, Bologna, Italy' and it should have been given as 'Istituto Scientifico Romagnolo per lo Studio e la Cura dei Tumori (IRST) IRCCS, Meldola, Italy and not Department of Experimental, Diagnostic and Specialty Medicine, University of Bologna, 40138, Bologna, Italy.'Furthermore, the original version of this Article contained an error in the spelling of the authors Alberto L'Abbate and Pietro D'Addabbo, an acute accent was used instead of an apostrophe for these authors names.These errors have now been corrected in both the PDF and HTML versions of the Article.

Frequent NRG1 fusions in Caucasian pulmonary mucinous adenocarcinoma predicted by Phospho-ErbB3 expression.

NRG1 fusions were recently reported as a new molecular feature of Invasive Mucinous Adenocarcinoma (IMA) of the lung. The NRG1 chimeric ligand acts as a strong inductor of phosphorylation and tyrosine kinase activity of the ErbB2/ErbB3 heterodimer, thus enhancing the PI3K-AKT/MAPK pathways. The NRG1 fusions were widely investigated in Asian IMA cohorts, whereas just anecdotal information are available about the occurrence of NRG1 fusions in IMA Caucasian population. Here we firstly explored a large Caucasian cohort of 51 IMAs and 34 non-IMA cases for the occurrence of NRG1 rearrangements by fluorescent in situ hybridization (FISH) and RNA target sequencing. FISH results were correlated to the immunohistochemical expression of phosphorylated-ErbB3 (pErbB3) receptor and the mutational status of KRAS, EGFR and ALK genes. The NRG1 rearrangements were detected in 31% IMAs and 3% non-IMAs and the CD74-NRG1 fusion transcript variant was characterized in 4 NRG1-positive IMAs. Moreover, pErbB3 expression was found to be strictly associated to the mucinous pattern (p = 0.012, Chi-square test) and all IMA cases showing aberrant expression of pErbB3 demonstrated NRG1 rearrangements. No significant correlation between NRG1 rearrangements and EGFR, KRAS or ALK mutations respectively, was observed. We report for the first time that NRG1 fusions are driver alterations clearly associated with mucinous lung adenocarcinoma subtype of Caucasian patients and not exclusive of Asiatic population. pErbB3 immunostaining may represent a strong predictor of NRG1 fusions, pointing out the detection of pErbB3 by IHC as a rapid and effective pre-screening method to select the NRG1-positive patients.

MYC-containing amplicons in acute myeloid leukemia: genomic structures, evolution, and transcriptional consequences.

Double minutes (dmin), homogeneously staining regions, and ring chromosomes are vehicles of gene amplification in cancer. The underlying mechanism leading to their formation as well as their structure and function in acute myeloid leukemia (AML) remain mysterious. We combined a range of high-resolution genomic methods to investigate the architecture and expression pattern of amplicons involving chromosome band 8q24 in 23 cases of AML (AML-amp). This revealed that different MYC-dmin architectures can coexist within the same leukemic cell population, indicating a step-wise evolution rather than a single event origin, such as through chromothripsis. This was supported also by the analysis of the chromothripsis criteria, that poorly matched the model in our samples. Furthermore, we found that dmin could evolve toward ring chromosomes stabilized by neocentromeres. Surprisingly, amplified genes (mainly PVT1) frequently participated in fusion transcripts lacking a corresponding DNA template. We also detected a significant overexpression of the circular RNA of PVT1 (circPVT1) in AML-amp cases versus AML with a normal karyotype. Our results show that 8q24 amplicons in AML are surprisingly plastic DNA structures with an unexpected association to novel fusion transcripts and circular RNAs.

The Hidden Genomic and Transcriptomic Plasticity of Giant Marker Chromosomes in Cancer.

Genome amplification in the form of rings or giant rod-shaped marker chromosomes (RGMs) is a common genetic alteration in soft tissue tumors. The mitotic stability of these structures is often rescued by perfectly functioning analphoid neocentromeres, which therefore significantly contribute to cancer progression. Here, we disentangled the genomic architecture of many neocentromeres stabilizing marker chromosomes in well-differentiated liposarcoma and lung sarcomatoid carcinoma samples. In cells carrying heavily rearranged RGMs, these structures were assembled as patchworks of multiple short amplified sequences, disclosing an extremely high level of complexity and definitely ruling out the existence of regions prone to neocentromere seeding. Moreover, by studying two well-differentiated liposarcoma samples derived from the onset and the recurrence of the same tumor, we documented an expansion of the neocentromeric domain that occurred during tumor progression, which reflects a strong selective pressure acting toward the improvement of the neocentromeric functionality in cancer. In lung sarcomatoid carcinoma cells we documented, extensive "centromere sliding" phenomena giving rise to multiple, closely mapping neocentromeric epialleles on separate coexisting markers occur, likely due to the instability of neocentromeres arising in cancer cells. Finally, by investigating the transcriptional activity of neocentromeres, we came across a burst of chimeric transcripts, both by extremely complex genomic rearrangements, and cis/trans-splicing events. Post-transcriptional editing events have been reported to expand and variegate the genetic repertoire of higher eukaryotes, so they might have a determining role in cancer. The increased incidence of fusion transcripts, might act as a driving force for the genomic amplification process, together with the increased transcription of oncogenes.

t(15;21) translocations leading to the concurrent downregulation of RUNX1 and its transcription factor partner genes SIN3A and TCF12 in myeloid disorders.

Through a combined approach integrating RNA-Seq, SNP-array, FISH and PCR techniques, we identified two novel t(15;21) translocations leading to the inactivation of RUNX1 and its partners SIN3A and TCF12. One is a complex t(15;21)(q24;q22), with both breakpoints mapped at the nucleotide level, joining RUNX1 to SIN3A and UBL7-AS1 in a patient with myelodysplasia. The other is a recurrent t(15;21)(q21;q22), juxtaposing RUNX1 and TCF12, with an opposite transcriptional orientation, in three myeloid leukemia cases. Since our transcriptome analysis indicated a significant number of differentially expressed genes associated with both translocations, we speculate an important pathogenetic role for these alterations involving RUNX1.

Therapeutic implications of intratumor heterogeneity for TP53 mutational status in Burkitt lymphoma.

Therapeutic implications of intra-tumor heterogeneity are still undefined. In this study we report a genetic and functional analysis aimed at defining the mechanisms of chemoresistance in a 43-year old woman affected by stage IVB Burkitt lymphoma with bulky abdominal masses and peritoneal effusion. The patient, despite a transient initial response to chemotherapy with reduction of the bulky masses, rapidly progressed and died of her disease. Targeted TP53 sequencing found that the bulky mass was wild-type whereas peritoneal fluid cells harbored a R282W mutation. Functional studies on TP53 mutant cells demonstrated an impaired p53-mediated response, resistance to ex vivo doxorubicin administration, overexpression of DNA damage response (DDR) activation markers and high sensitivity to pharmacologic DDR inhibition. These findings suggest that intra-tumor heterogeneity for TP53 mutational status may occur in MYC-driven cancers, and that DDR inhibitors could be effective in targeting hidden TP53 mutant clones in tumors characterized by genomic instability and prone to intra-tumor heterogeneity.

Genomic organization and evolution of double minutes/homogeneously staining regions with MYC amplification in human cancer.

The mechanism for generating double minutes chromosomes (dmin) and homogeneously staining regions (hsr) in cancer is still poorly understood. Through an integrated approach combining next-generation sequencing, single nucleotide polymorphism array, fluorescent in situ hybridization and polymerase chain reaction-based techniques, we inferred the fine structure of MYC-containing dmin/hsr amplicons harboring sequences from several different chromosomes in seven tumor cell lines, and characterized an unprecedented number of hsr insertion sites. Local chromosome shattering involving a single-step catastrophic event (chromothripsis) was recently proposed to explain clustered chromosomal rearrangements and genomic amplifications in cancer. Our bioinformatics analyses based on the listed criteria to define chromothripsis led us to exclude it as the driving force underlying amplicon genesis in our samples. Instead, the finding of coexisting heterogeneous amplicons, differing in their complexity and chromosome content, in cell lines derived from the same tumor indicated the occurrence of a multi-step evolutionary process in the genesis of dmin/hsr. Our integrated approach allowed us to gather a complete view of the complex chromosome rearrangements occurring within MYC amplicons, suggesting that more than one model may be invoked to explain the origin of dmin/hsr in cancer. Finally, we identified PVT1 as a target of fusion events, confirming its role as breakpoint hotspot in MYC amplification.

FOXP1 and TP63 involvement in the progression of myelodysplastic syndrome with 5q- and additional cytogenetic abnormalities.

BACKGROUND: The progression of low-risk del(5q) myelodysplastic syndrome to acute myeloid leukemia is increased when associated with mutations of TP53, or with additional chromosomal abnormalities. However, to date the prognostic impact and molecular consequences of these rearrangements were poorly investigated. Single additional alterations to del(5q) by balanced chromosome rearrangements were rarely found in myelodysplasia. In particular, balanced alterations involving TP63 and FOXP1 genes were never reported in the literature. CASE PRESENTATION: Here we report on a 79-year woman with an aggressive form of myelodysplastic syndrome with del(5q), no TP53 mutation, and a novel complex rearrangement of chromosome 3 in bone marrow cells. Our results revealed that the FOXP1 and TP63 genes were both relocated along chromosome 3. Strikingly, immunohistochemistry analysis showed altered protein levels, disclosing that this rearrangement triggered the expression of FOXP1 and TP63 genes. FOXP1 was also found activated in other patients with myelodysplasia and acute myeloid leukemia, showing that it is an important, recurrent event. CONCLUSIONS: We document an apparent role of FOXP1 and TP63, up to now poorly documented, in the progression of MDS in our patient who is lacking mutations in the TP53 tumor suppressor gene normally associated with poor outcome in myelodysplastic syndrome with 5q-. Finally, our results may suggest a possible broader role of FOXP1 in the pathogenesis and progression of myelodysplasia and acute myeloid leukemia.

Proteolysis of MOB1 by the ubiquitin ligase praja2 attenuates Hippo signalling and supports glioblastoma growth.

Human glioblastoma is the most frequent and aggressive form of brain tumour in the adult population. Proteolytic turnover of tumour suppressors by the ubiquitin-proteasome system is a mechanism that tumour cells can adopt to sustain their growth and invasiveness. However, the identity of ubiquitin-proteasome targets and regulators in glioblastoma are still unknown. Here we report that the RING ligase praja2 ubiquitylates and degrades Mob, a core component of NDR/LATS kinase and a positive regulator of the tumour-suppressor Hippo cascade. Degradation of Mob through the ubiquitin-proteasome system attenuates the Hippo cascade and sustains glioblastoma growth in vivo. Accordingly, accumulation of praja2 during the transition from low- to high-grade glioma is associated with significant downregulation of the Hippo pathway. These findings identify praja2 as a novel upstream regulator of the Hippo cascade, linking the ubiquitin proteasome system to deregulated glioblastoma growth.

Thrombocytopenia-absent-radius syndrome in a child showing a larger 1q21.1 deletion than the one in his healthy mother, and a significant downregulation of the commonly deleted genes.

Thrombocytopenia-absent-radius (TAR) syndrome is a rare condition characterized by thrombocytopenia and bilateral absence of the radii with presence of both thumbs. The phenotype has a variable expression. A 200 kb minimal deletion at 1q21.1 is present in all patients. However, the microdeletion, ranging up to 1100 kb in length, is not sufficient to cause the disease. Indeed it is present in 75-80% of unaffected parents. It is assumed that the phenotype develops only in the presence of one or more additional, as-yet-unknown, deletion modifiers (mTARs). We report here on a child affected by TAR syndrome associated with Langerhans cell histiocytosis. Unexpectedly, he showed a 2.029 kb deletion at 1q21.1, almost twice that of the unaffected mother (957 kb). Interestingly, the mother-to-son increased size of the deleted region was already observed in two cases of constitutional diseases, although both resulting as chromosomal terminal deletions. Noteworthy, qPCR experiments, never before performed for patients with TAR syndrome, disclosed that the proband had a statistically significant downregulation of the majority of the genes mapping inside the part of the deletion shared with the mother. The mother, on the contrary, did not show the same downregulation. In summary, the present report adds new insights on the pathogenesis of TAR syndrome, that may represent fruitful directions for future research.